Excitation spectrum of d-wave Fermi surface deformation 
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Several instabilities competing with the d-wave singlet pairing were proposed for high-T c cuprates. 
One of them is the d- wave Fermi surface deformation (dFSD), which is generated by forward scat- 
tering. In this paper, correlation functions of the dFSD are calculated within the random phase 
approximation. In the normal state, the excitation spectrum shows a low energy peak, which 
smoothly connects to critical fluctuations of the dFSD at lower temperature. The competition with 
the d-wave pairing, however, blocks the critical fluctuations. The whole spectral weight is trans- 
ferred to high energy and a pronounced peak appears there in the d-wave pairing state. This peak 
is an overdamped collective mode of the dFSD and can grow to be a resonance mode at moderate 
finite wavevectors. 

PACS numbers: 71.18.+y, 71.10.Fd, 74.25.-q, 74.72.-h 
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High-T c cuprates are doped Mott insulators. The 
parent compounds are antiferromagnetic Mott insula- 
tors, which become high-T c superconductors with car- 
rier doping. The superconducting state does not have an 
isotropic gap such as BCS superconductors, but has an 
anisotropic gap with the <i-wave symmetry. 

It has been recognized that there may be several in- 
stabilities competing with the d-wave superconductiv- 
ity. Effects of the antiferromagnetism may play a cru- 
cial role still in the d-wave superconducting state, and 
their competition was discussed using a concept of the 
SO(5) symmetry. 0,0 Another idea, a self-organized one- 
dimensional charge order in the CuC>2 plane — spin- 
charge stripes hypothesis, was proposed to discuss several 
experimental data[3|, |^ and phenomenological theories 
were developed. [j| Through a microscopic analysis of the 
two-dimensional (2D) t-J model by 1/N expansions for 
Hubbard operators, the d-wave charge density order was 
proposed. (3 This phase has bond currents forming stag- 
gered flux, and was discussed in different contexts. 0, 
A similar state, called the staggered flux phase, was 
proposed in the SU(2) slave-boson formalism in the t-J 
model. (3 In this scheme, the exact SU(2) gauge symme- 
try at half-nlling[ljj was invoked also at finite doping; the 
underlying theoretical concept is quite different from the 
d-wave charge density order. All these possible compet- 
ing orders come from electron-electron correlations with 
large momentum transfer near q = (vr, 7r). 

Recently another competing order was proposed, 
ll2T | the d-wave Fermi surface deformation (dFSD).(l^ 
The Fermi surface (FS) expands along the fc^-direction 
and shrinks along the /c^-direction (or vice versa). This 
order has to be distinguished from the above possible 
competing orders. The channel of this instability is for- 
ward scattering with q = (0,0). The dFSD was first dis- 
cussed for the 2D t-J model^lj and Hubbard model. (l2^ 
It was tested in several renormalization group schemes 
applied to the Hubbard model. [H Q The dFSD 
was investigated also in perturbation theories for the 
Hubbard model 0, [HI, in the mean-field theory for the 



extended Hubbard model, [lg and in phenomenological 
models. 0, |2(], In the continuum (not lattice) model 
FS deformation was investigated in analogy to the ne- 
matic phase in liquid crystals. p2ll23| 

In accord with a result in the Hubbard model, the 
analysis of the t-J modeljllj showed that an instability 
of the dFSD competed with a more dominant instabil- 
ity, the d-wave singlet pairing, and was usually masked 
and not seen. However, it was shown that the presence 
of a small extrinsic anisotropy was sufficient to mani- 
fest the dFSD. This implies that while the sponta- 
neous instability of the dFSD does not take place, the 
electron system still has an appreciable susceptibility of 
the dFSD and is sensitive to the external anisotropy; 
the FS is softened. This idea was invoked for LSCO 
systems through the consideration of band parameter 
dependences^]] and magnetic excitation spectra. |24| 

In this letter, we investigate dynamical properties of 
the dFSD. Since the instability of the dFSD is signaled 
by divergence of its static susceptibility at q = 0, we 
focus on the dynamical susceptibility near q = and cal- 
culate it within the random phase approximation (RPA). 
In the normal state the excitation spectrum shows a low 
energy peak, which smoothly connects to critical fluc- 
tuations of the dFSD at lower temperature. The criti- 
cal fluctuations are, however, blocked by the more domi- 
nant d-wave pairing instability. The low energy spectral 
weight is suppressed and vanishes at zero temperature. 
Instead the spectral weight is transfered to high energy 
and we find a pronounced peak there. This peak is an 
overdamped collective mode of the dFSD and can grow 
to be a resonance mode at moderate finite wavevectors. 

To investigate correlations of the dFSD, we take the 
2D t-J model on the square lattice, 
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defined in the Fock space with no doubly occupied sites. 
Here Ci a (Si) is an electron (a spin) operator. The 
is the Ith (I < 2) neighbor hopping integral, and we de- 
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note = t, i (2) = t'. The J(> 0) is the superex- 
change coupling between nearest-neighbor sites. We in- 
troduce U(l) slave-particles as c\ a = fj a bi, where /$„ 
(bi) is a fermion (boson) operator that carries spin a 
(charge e), and Si = \ f} a er a p fi p with Pauli matrix 
er. This is an exact transformation. We then decou- 
ple the interactions with the so-called resonating- valence- 
bond mean fields: Xt=(2 ct fLfi+ra), {b\b i+T ), and 
A T =(fiifi+Ti - fiifi+rf}, where r = r i - n denotes 
the direction. These mean fields are assumed to be real 
constants independent of sites i. We approximate the 
boson to condense at the bottom of its band, which is 
reasonable at low temperature T, and obtain the follow- 
ing Hamiltonian: 
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with a global constraint J2a(ficrfi>?) = 1 — ^- Here = 
—2 (F x cos k x + F y cos k y + 2t'5 cos k x cos k y ) — ^, — 
-\j {A x cosk x + A v cos ky), and F x{y ) = tS + \JXx{y), 
with 8 being hole density and /j, the chemical potential. 
The mean fields are determined self-consistently by mini- 
mizing the free energy. The isotropic state \ x — \y is sta- 
bilized and the ci- wave singlet pairing, A x = —A y = A , 
sets in at low T. 

The advantages of this formalism[2£j are (i) the phase 
diagram on T versus S catches essential physics of high-T c 
cuprates, (ii) magnetic excitation in actual systems was 
consistently described on the basis of fermiology includ- 
ing material dependence.pl . IM |27| and (iii) the dFSD 
channel was shown to exist in the J-term.|ll| which en- 
able us study its competition with the (i-wave pairing on 
an equal footing. 

To analyze correlations of the dFSD, we define d-wave 
weighted fermion density, 
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with dk — ^(cosfca; — cos k y ). The spontaneous dFSD is 
described by (x<j(0)) 7^ 0, which is, however, prohibited 
by competition with the d-wave singlet pairing. In the 
d-wave singlet state, fluctuations of xAq) induce fluc- 
tuations of extended s-wave pairing. 11] To include this 
effect wc also define 



A s (q) = X>fc(/fe+!T/-fe+fl-/fe-! T /-fe-§ J ' ( 4 ) 



with Sk = ^(cosfca, + cos ky). Thus the correlation func- 
tion forms a 2 x 2 matrix, 
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where n^{q,uj) = Kq 1 (<7,w), and 



l ^ 1 {q,w) = -J dte^+^iiUq^^M-q)}}*), (6) 
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K J 2 (g, w) = - / dte^ +ir »(lMq, t),A s (-q)}) , (7) 
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dte^ +ir ^([A s (q, t),A s (~q)}) . (8) 



Here x d (q, t) = c lH " 1 Xd(q)e- iHot , and A s (q, t) = 
e iH t A s (q)e~~ lHot . The bracket (• ■ • }o denotes an expec- 
tation value under the Hamiltonian J2J), and [• , •] is the 
commutator; N is the total number of lattice sites. We 
consider interactions in the RPA, 
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In this letter, we focus on K n (q, uj) with q w 0, and 
investigate its spectral weight in both the normal state 
and the d-wave singlet pairing state. Full results includ- 
ing other components of k will be shown elsewhere, psj 
We take band parameters, t/J — 4 and t' jt = —1/6, for 
which the dFSD is known to be prominent : |l l| a doping 
rate is fixed to S — 0.10. In Eqs.©-©, the value of T is 
a positive infinitesimal and we take T = 10~ 4 J (Fig. ^) 
or r = 0.01 J (Fig. [21 in numerical calculations. 

Figure QJa) shows lm.K 11 (q, uj) as a function of uj 
for several choices of T in the normal state at q = 
(0.01 x 2-7T, 0). The spectral weight concentrates at 
low energy for all T. This is due to a property of 
ImKQ 1 (q, uj) that particle-hole excitations obey the re- 
lation, uj = £fc +q /2 — Cfc-q/2! w becomes small for small 
q. The low energy spectral weight increases with de- 
creasing T and the peak position shifts closer to zero 
energy. This enhancement comes from the interactions 
in the RPA, Eq. 10, and smoothly connects to critical 
fluctuations of the dFSD at lower T. While the present 
lowest temperature (T = 0.13 J) is much higher than the 
critical temperature of the dFSD (T^fsd = 0.038 J), the 
low energy weight of Im/i n (g, uj) substantially increases 
in comparison with Im/tJ 1 (q, uj) as shown in the inset 
of Fig. nja). In this sense, the enhancement of the low 
energy peak is a precursor of collective fluctuations of 
the dFSD. In Fig. ^b), we plot ImK 11 (q, uj) for several 
choices of q (|| [10]) at given T. The spectral weight 
spreads to higher energy with increasing |q|, but the peak 
position stays at relatively low energy, which is due to the 
enhancement by the RPA. In Fig.^c), we summarize the 
excitation spectrum on the plane of uj vs. q. The shaded 
region is a gapless particle-hole continuum and the up- 
per edge increases linearly with q. The peak energy of 
IuiK 11 (q 1 uj) disperses linearly with q at low q within nu- 
merical accuracy. The gradient of the g-linear becomes 
small at low T, which is due to the enhancement of low 
energy fluctuations of the tfFSD. These qualitative fea- 
tures have been checked also along q \\ [11]. 
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FIG. 1: Normal state, (a) uj dependence of IniK 11 (q, uj) for 
several choices of T at q — (0.01 x 2ir, 0). The result at 
T = 0.13J is compared with Itcikq (q, uj) in the inset, (b) 
uj dependence of ImK 11 (q, u>) for several choices of q along 
the [10] direction, (c) Excitation spectrum on the plane of 
u> vs. q. The shaded region is a particle-hole continuum. 
The open circle (triangle) corresponds to the peak energy of 
ImAt 11 (q, uj) at a given q at T = 0.15 J (0.30 J); q is defined 
as q = (q x 2iv, 0). 



FIG. 2: d-wave pairing state, (a) uj dependence of S 11 (q, uj) 
for several values of T at q = (0.01 x 2ir, 0). (b) uj dependence 
of ImK 11 (q, uj) for several choices of q along the [10] direction 
(T = 0.01J). (c) Excitation spectrum on the plane of uj vs. 
q. The shaded region is a continuum of excitations. The open 
circle corresponds to the peak energy of ImK 11 (q, uj), uj TC b, at 
a given q (T — 0.01J); q is defined as q — (q x 2-k, 0). i2) re s is 
estimation by Eq. I1UI . 



Further decreasing T below values shown in Fig.^a), 
the d-wave singlet pairing instability takes place, which 
competes with the g?FSD and prohibits the spontaneous 
dFSD. This competition is shown in Fig.[5Ja); we take the 
spectral function, S 11 (q, uj) — 2ImK 11 (q, w)/(l — e~ w / T ), 
to see low energy structures also on the same scale, and 
plot its uj dependence for several choices of T. The low 
energy weight is suppressed with decreasing T and is 
transferred to higher energy to form a second peak there 
(see the results for T > 0.12 J). While the low energy 
weight vanishes at T = 0, the second peak grows to be a 
pronounced peak at low T. To see its dispersive features, 
we calculate the uj dependence of ImK 11 (tj, uj) for several 
choices of q at low T. Figure Hfb) shows that the peak 
width becomes narrower with |q| and a sharp peak ap- 
pears at moderate |q| (> 0.05 x 27r). This is a resonance 
peak and an in-gap collective mode of the dFSD, namely 
a bound state (the finite peak width of the resonance 
is due to r > in the numerical calculations). To see 
this, we calculate the gap energy of Im/tJ (q, uj), namely 
a lower edge of a continuum of excitations, and plot it as 
a function of q as well as the peak energy of luiK 11 (q, uj), 
uj ics , in Fig. H[c). The lower edge increases with q and 
u) Tes is located inside the gap (in q > 0.05), which gives 
rise to the resonance peak. This resonance does not ap- 
pear along q \\ [11] at least up to q = (0.10x27r, 0.10x27r) 



because of the extention of the continuum spectrum down 
to lower energy. 

To understand the dispersion relation of the reso- 
nance, we first approximate n ll {q, uj) given by Eq. @ 
to K n (q, uj) = n^iq, uj)/(1 — 3JnQ 1 (q, uj)/2); we have 
checked numerically that quantitative changes by this ap- 
proximation are not noticeable. The dispersion is then 
determined by the condition, 

l-pR CK 1 1 {q,uj ICS ) = 0. (10) 

This equation can have two solutions for a given q and 
the uj res is the smaller one. We solve Eq. (|10l) numerically 
and compare cl> res with uj les in Fig. EKc). We see a good 
agreement in a moderate g-region, where the resonance 
appears. A poor agreement in a small (7-region is due 
to finite weight of lm.KQ 1 (q, uj), which invalidates using 
Eq. (jlOJl . We, however, see that the dispersive features 
are well characterized by Eq. I|10|) in a whole g-region in 
Fig. Etc). 

It should be noted that Eq. (|10fl has solutions for any 
q shown in Fig.[2tc). This is due to the enhancement of 
ReKQ 1 (t7, uj) by the d-wave form factor in Eq. Q. Since 
Eq. QlOjl describes an in-gap collective mode at moder- 
ate q, the peak of ImK 11 (t7, uj) in the small (7-region is 
regarded as an overdamped collective mode of the cfFSD. 

We have investigated the RPA excitation spectrum of 
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the <iFSD within the slave-boson mean-field approxima- 
tion to the 2D t-J model. In a normal state, excitation 
spectrum shows a low energy peak, which connects to 
critical fluctuations of the dFSD at lower T. The compe- 
tition with the ci-wave pairing, however, blocks the criti- 
cal fluctuations. The whole spectral weight is transferred 
to high energy and a pronounced peak appears there in 
the rf-wave pairing state. This peak is an overdamped 
collective mode of the g?FSD and can grow to be a reso- 
nance mode at moderate q. 

While these results are obtained in the RPA, we expect 
that higher order corrections will not modify appreciably 
at least the results near T = (Fig.[3J), since the boson 
condenses at the bottom of its band and the U(l) gauge 
field describing fluctuations around the mean fields is not 
relevant. On the other hand, the results of Fig. are 
obtained at finite T and a g-linear behavior of the low 
energy peak might not be a robust property. 

Correlations of the dFSD are ingredients of both the 
2D t-J modeldJ and Hubbard model[H Q E E3. 
Their implications for high-T c cuprates are interesting. 
Since appreciable correlations of the g?FSD make the 
electron system sensitive to an extrinsic anisotropy be- 
tween the x-direction and the y-direction, even a small 
anisotropy can be sufficient to lead to the dFSD, possi- 
bly a quasi-lD FS in each CuC>2 plane. This possibility 
was proposed for Nd-doped LSCO systems. [ll| Fluctua- 
tions of the dFSD in such a quasi-lD state will be investi- 
gated elsewhere. 28] In the absence of a (static) xy-spatial 
anisotropy, the present theory is applicable and we expect 
the excitation spectrum Fig.^c) in the normal state and 
Fig- Etc) in the d-wave pairing state.[2^ As a direct test, 



however, conventional optical methods are not sufficient, 
since they measure a quantity with q = and will not 
reach a finite (/-region, especially the region where the 
resonance mode is predicted. Indirectly, searching some 
phonon anomalies may be promising since fluctuations of 
the <£FSD are expected to couple with a lattice degree of 
freedom. However, we do not have calculations on such 
coupled systems at present. 

Fluctuations of the dFSD should not be confused with 
those of spin-charge stripes. 0, 13 W The dFSD is gen- 
erated by forward scattering while formation of spin- 
charge stripes requires an interaction with large momen- 
tum transfer. The underlying physics is different, (ii) 
Fluctuations of the cZFSD are relevant in systems near the 
breaking of the square lattice symmetry while stripe fluc- 
tuations make sense in systems near translational sym- 
metry breaking. 

While some hidden orders are often discussed in the 
connection with the pseudogap.|30l| correlations of the 
<iFSD are not related directly to the pseudogap in the 
slave-boson scheme. [2j| However, their effects may con- 
tribute to pseudogap behaviors additively, since the 
present FS deformation has the ci-wave symmetry, the 
same symmetry as the pseudogap. 
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